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^ ■ ABSTRACT 

m 

The daily variation of the solar photocenter over some 11 years is derived from 

Pm ! the Mount Wilson data reprocessed by Ulrich et al. 2010 to closely match the 

. ' surface distribution of solar irradiance. The standard deviations of astrometric 

p ' jitter are 0.52 /lAU and 0.39 /lAU in the equatorial and the axial dimensions, 

respectively. The overall dispersion is strongly correlated with the solar cycle, 

reaching 0.91 /lAU at the maximum activity in 2000. The largest short-term 



o3 \ deviations from the running average (up to 2.6 /xAU) occur when a group of 
large spots happen to he on one side with respect to the center of the disk. The 

^ \ amplitude spectrum of the photocenter variations never exceeds 0.033 //AU for 

00 \ the range of periods 0.6-1.4 yr, corresponding to the orbital periods of planets 

rr I in the habitable zone. Astrometric detection of Earth-like planets around stars 

'sf \ as quiet as the Sun is not affected by star spot noise, but the prospects for more 

in 

o 

o 



•n I active stars may be limited to giant planets. 



Subject headings: open clusters and associations: general — stars: kinematics 
binaries: general 



1. Introduction 

The prospects of finding habitable planets orbiting nearby solar-type stars are to a large 
degree associated with the ultra-precis e astrometric instruments under deyelopment or con- 



struc tion, such a s the SIM Observatory (IShao et al.ll2009l : lUnwin et al.ll2008t ICatanzarite et al 



20061 ) and Gaia (jCasertano et al.ll2008l ). The Earth orbiting the Sun produces an observable 
astrometric wobble of 3 /xAU (micro-AU) and a radial velocity variation of 0.089 m s~^, if 
seen equator-on (inclination i = 90°). A pole-on configuration {i = 0° or 180°) is optimal for 
astrometry, because the reflex motion signal is present in both dimensions of the sky projec- 
tion, whereas the radial velocity amplitude, which is proportional to sini, drops to zero. For 
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a high signal-to-noise (S/N) ratio detection, the total error budget of the prospective exo- 
planet detection techniques should be well below these values. The astrophysical jitter caused 
by the rotation of magnetic features on the surface (spots, faculae) potentially becomes a 
significant part of the total observational error in this domain of accuracy, even for fairly 
common inactive dwarfs. The recently completed double blind test of planet-detection capa- 
bilities w ith the SIM Obse rvatory produced encouraging results even for complex multiple 
systems (JTraub et al.l 120091 ). but magnetic jitter was not taken into account in those simu- 
lations. Semi-analytical considerations in ( JMakarov et al.l 120091 ) supported by Monte-Carlo 
simulations and indirect observational evidence determined that the astrometric method is 
much less sensitive to the effects of magnetic features, at least by an order of magnitude 
for solar-type stars, tha n the Doppler shift technique. Similar conclusions were drawn by 
Catanzarite et al.l (|2008[ ) for the Sun, based on a sophisticated model of sunspot activity and 
extensive numerical simulation. 

The variable distribution of surface brightness is the cause of both the astrometric jitter 
and the total flux variation. Therefore, the former can be estimated from the latter, given 
a model describing the number of spots (or bright features), their size and lifetimes. The 
dispersion of the total solar irra diance (TSI) observed for a few decades with a number of 
sateUites (JFrohlich &: LeanI Il998l ) sets the standard of quiet stars (4 ■ 10~^, or 400 ppm in 
normalized flux). The light curves from Kepler for some 120000 stars rev ealed that at leas t 
half of all solar- type dwarfs in the field are at least as quiet as the Sun ( iBasri et al.ll2010l ). 
displaying lower levels of photometric dispersion. By extrapolation, about half of nearby 
dwarfs should be as amenable to exoplanet detection as the Sun. The aim of this paper is 
to determine the solar astrometric jitter directly from observations. The results can be used 
to verify and correct the existing models of spot activity for other stars. 



The data 



We made use of the set of 2881 images of the solar disk reconstructed from Mount 
Wilson Observator y magnetogram and intensity ratio images for the period 1996 to 2007 
( lUlrich et al.l l2010l ). As described in that paper, two options of the AutoClass analysis 
yielded an 18 and a 37 class set of solar surface features. The original Mount Wilson images 
in the intensity ratio /a525o/-^A5237 and magnetic field strength l-Bx^^.^ nl were expanded in th e 
corresponding sets of classes and calibrated against the TSI data (JFrohlich fc LeanI Il998l ). 
The resulting solar irradiance images have the full resolution of the original data, with a value 
of normalized TSI assigned to each pixel. Essentially, these images are two-dimensional maps 
of bolometric surface intensity. The reconstructed images are only 95% of the full solar disk 
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in radius. The units of pixel values are W-m ^, and they are normalized to unit total area. 
For this study, we used the 37-class reconstruction, which provides a sl ightly better match 



to the TSI, although this improvement is not important for our purpose. lUlrich et al.l (120 lOf ) 
discuss the "ring effect" in the reconstructed TSI maps, which manifests itself as a kind 
of diffraction ring surrounding high-contrast, compact features, especially dark spots. The 
spatial scale of these artefacts is small, and the impact on the estimated magnetic jitter is 
negligible. 

The processing of the TSI maps was straightforward. Each of the 2881 images was 
integrated in first moment with respect to a fixed pixel row, e.g.. 



LAX / ^ J {•^ijj Vij ) -^ij/ / J -' \-^iji Vijii 



and similarly for Ay. The x axis is aligned with the solar equator, and the y axis with 
the rotation axis. The derived offsets were brought to zero mean and converted to /lAU. 
Since the TSI images are based on the observed distributions of magnetic strength and the 
ratio of two monochromatic intensities, they are almost devoid of the limb darkening effects. 
Therefore, we do not apply any limb darkening corrections as the counter-acting effects of 
the broad-band intensity reduction to ward the limb and the increasing contrast of faculae 



with respect to the local photosphere (JFoukal et al.l 120041 ) have been balanced in the images 
by construction. 



3. Astrometric jitter from the magnetic activity 

The main results of this study are shown in Fig. [H The top panel displaying the 
integrated TSI is shown only for reference, because the TSI images were constructed to match 
the space-based TSI measurements as close as possible, thus, the plot does not contain any 
new information. The middle and bottom panels show the daily photocenter offsets of the 
solar disk over some 11 years in the equatorial and axial dimensions, respect ively. The TSI 



has b een for some time known to vary systematically with the solar cycle (jWillson et al 



198ll ). but the exact contribution of the opposing effects (faculae and spots) to this process 
is still somewhat controversial. The Sun is brighter when its level of magnetic activity 
i s elevated, while y ounger and more active stars tend to display the opposite correlation 



(JRadick et al.l Il998l ). The peak-to-peak variation of TSI is approximately 2 W-m~^, or 
0.15%. The maximum of activity in 2000-2002 is marked with a brighter average and a 
much larger dispersion of TSI. A running average standard deviation over 3 months varies 
between 0.14 yuAU and 0.91 /lAU in Ax. The overall standard deviations are 0.52 yuAU for 
Ax and 0.39 yuAU for Ay. Predictably, the dispersion of Ay is smaller, because major spots 
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occur in the equatorial zone of the disk. 
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Fig. 1. — Daily variations of the total solar irradiance, and the photocenter position in the 
equatorial (Ax) and axial (Ay) dimensions, derived from the reconstructed TSI maps. 

The distinct variation in Ay with a period of 1 year is probably an artefact related to 
the seasonal change of the solar obliquity. In September, the solar axis is tilted by 7.25°, with 
its north pole toward the Earth. Since the distribution of magnetic features is symmetric 
about the solar equator, more spots and faculae appear in the southern part of the disk than 
in the northern part, as seen by a ground-based observer. This 1-year modulation of the Ay 
offset is germane to the Sun and will not be observed on other stars. A periodogram analysis 
determined that the best fitting sinusoid has a period of 1 yr and an amplitude of 0.11 /lAU. 

The running average of the Ax perturbations is not correlated with the solar cycle, 
because the probability of magnetic features to be located in the eastern or western hemi- 
spheres is the same. The scatter of daily variation is, however, strongly correlated with the 
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level of magnetic activity. The distribution of TSI perturbations around the running mean 
is asymmetric, with a longer tail extending toward fainter intensities. These short-lived dips 
in TSI are caused by particularly large sunspot groups, and normally happen at times of 
elevated activity. The distribution of astrometric perturbations around the running mean 
are fairly symmetric, as expected. The largest astrometric perturbation Ax took place on 
2000.09.19, when a distant observer would have detected a westward photocenter shift of 
2.6 /iAU. Fig. |2] depicts the distribution of surface brightness of the Sun on that date. The 
color scale is arbitrary, adjusted in such a way that the dark spots (black and blue) and the 
faculae (yellow) can be clearly seen. This unusually large perturbation was caused by two 
sunspot groups, which happened to be in the eastern hemisphere, and a bright active region 
in the opposite western hemisphere. 




Fig. 2. — Solar surface irradiance distribution on 2000.09.19, color-coded to emphasize the 
dark spots and the bright plage areas. North is up, east is left. An exceptionally large 
deviation of —2.6 yuAU along the equator would be observed astrometrically. 



The temporal variations of solar surface brightness are commonly considered to be 
caused by small- and medium-scale magnetic features, such as photospheric spots and plages. 
Is it possible that some other, non-magnetic structures are present in the surface brightness 
distribution, which would not be captured b y the Mount Wilson magnetograms that served as 
the original data for the TSI reconstruction? lUlrich et al.l ( 120101 ) found that the reconstructed 
irradiance followed the actual TSI variations, observed with Virgo, to a correlation coefficient 
of 0.9625. This corresponds to a standard deviation of 28% relative to the variation of TSI, 
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if the error introduced by the proxy reconstruction is uncorrelated with the true hght curve. 
Although it is not possible to tell how this error is distributed across the solar disk, we 
condend that the true astrometric jitter can not be larger than our estimates by more than 
this amount, which is 4% when added in quadrature. The error is carried by many classes 
of different spatial scale, but only large scale perturbations count in the integration in first 
moment. Most of this error is probably confined to small-scale structures, which should 
cancel in integration. Known imperfections, such as the ringing features centered on the 
centers of activity, caused by too many plage classes, do not result in an appreciable image 
shift. Besides, some of the TSI variation not accounted for by the reconstruction is due to 
magnetic features appearing in the outer rim of the solar disk (5% of the radius), which is 
missing in our data. 



4. Implications for Earth-like planet detection 

Planets revolving around their host stars produce sinusoidal variations in the astromet- 
ric position. The main harmonic of these variations has the orbital period of the planet, 
and an amplitude defined by the mass of the star, the mass of the planet, and the orbital 
semimajor axis. Planets in eccentric orbits also produce higher-order harmonics (overtones) 
of smaller amplitudes. The spectroscopic method of exoplanet detection, based on precision 
measurements of st ellar radial velocity, utilizes the Lomb-Scargle periodogram analysis (e.g.. 



Fischer et al.ll2008l ). which is aimed at detectin g statistically significant sinusoidal variations 



in irregularly sampled data with a zero mean ( IScargld 119821 ). A s imilar "joint" pow'er spec - 



trum periodogram was suggested for 2D astrometric detection by ICatanzarite et al.l ( l2008l ). 
In order to evaluate the impact of magnetic jitter on detection of low-mass planets, we em- 
ploy in this paper a generalized amplitude spectrum analysis, which, unlike the Lomb-Scargle 
periodogram, includes the constant term in the set of fitting functions. A constant offset in 
position is a physical astrometric parameter, which can not be simply subtracted from the 
raw data. 

The linear problem 

[1, cos(ajt), sin(cjt)] s = Ax (2) 

is solved for a grid of frequencies u = 2tt/p, where t is the column vector of observation 
times, and Ax is the column vector of x-coordinate offsets. If s is the Least-Squares solution 
(a 3- vector), the generalized amplitude spectrum as a function of period is 



D{p) = ^JSl + §l (3) 

The interpretation of this function is straightforward: it is the amplitude of the best-fitting 
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sinusoid with the given period p and a free phase. A planet with an orbital period p can 
be confidently detected if the amplitude of its main harmonic is much larger than D{p). It 
corresponds to the square root of the "periodogram power" of the Lomb-Scargle periodogram 
commonly used in exoplanet detection. 

Fig. [3] shows the D{p) spectrum of the Ax series derived from the TSI maps. The 
calculation was restricted to the range of habitable zone around the Sun, which corresponds 
to, somewhat generously, 0.6-1.4 years in orbital period. The spectrum is distributed non- 
uniformly, so that white noise would not be a good model for this kind of perturbation. 
The largest sinusoidal variation has a period of ~ 1.1 yr and an amplitude of 0.033 fiAU. 
The Earth generates an astrometric wobble of the Sun of 3 fiAU, whose power is almost 
entirely in the main harmonic (1 year) due to the very low eccentricity. Thus, the Earth can 
be detected by astrometric means at a very comfortable S/N of greater than 90, as far as 
magnetic jitter is concerned, if a similar observational cadence is achieved. For a solar-like 
star at 10 pc, the amplitude spectrum of magnetic jitter is not greater than 0.0033 /ias, 
which is vanishingly small compared to the expected single measurement precision of SIM 
(> 1 /ias) or Gaia (> 8 yuas). 



5. Discussion 

We determined that the solar jitter caused by magnetic features is small and should 
not preclude astrometrists from detecting habitable planets as small as the Earth. A few 
additional considerations should be made when extrapolating this result to other stars or to 
actual astrometric observations. 

The ultra-precise photometric measurements wit h Kepler confirme d that about half of 



all solar-type dwarfs are as quiet as the Sun, or more (jBasri et al.ll2010[ ). The stars that are 
more active than the Sun have much larger filling factors, and their magnetic features may be 
long-lived (up to a few weeks). Both photomet ric and astrometric y ariability is proportional 



to the average total area Ag occupied by spots (JMakarov et al.ll2009l ). In the model currently 



adopted for planet detection simulations with SIM (to be published elsewhere), the total 
area is related to the index of chromospheric activity log -R^k through 

log A, = 12.468 + 1.75 log i?ijK, (4) 

which yields 6137 micro-stellar-hemispheres (MSH) for the Sun. As roughly half of nearby 



field stars are more chromospherically active than the Sun (iGray et al.ll2003l ). they are ex- 
pected to have larger spot groups and, therefore, jitter amplitudes. At logi^^K = —4.6, the 
expected area and the jitter are 5 times the solar values, and at log-R^j^ = —4.2, 25 times. 



0.04 



0.035 




0.9 1 1.1 

Period years 



1.4 



Fig. 3. — Amplitude spectrum of solar astrometric jitter in the equatorial dimension for the 
range of periods corresponding to the habitable zone. 
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This mo del gives fairly accur ate predictions for the active star n^ Ceti, carefully studied with 



MOST ( iWalker et al. 1120071 ). Furthermore, an additional factor of 2-6 increase may ensue 
from the limited duration of an astrometric mission. Both SIM and Gaia will nominally 
operate for 5 years. This is considerably shorter than the solar cycle. If a solar analog is 
observed at the height of its activity, the level of magnetic jitter may be higher than the 
multi-year average. Thus, magnetic jitter can indeed pose the natural limit to the sensitivity 
of the astrometric method at stars significantly more active than the Sun. 

The amplitude spectrum periodogram shown in Fig. [3] is derived from a high-frequency 
observational cadence of 2881 data points over 11 years. We expect only about 200 2D 
measurements to be taken for each target star with SIM Lite during its 5-year nominal 
mission. If the astrometric jitter is completely random and uncorrelated, the variance of the 
amplitude of a random harmonic in the spectrum is inversely proportional to the number of 
data points. In that case, the amplitude spectrum will rise by a factor of a/2881/200 = 3.8 
for SIM. On the other hand, if the magnetic activity underlying the astrometric jitter is 
systematic and deterministic, the corresponding parts of the spectrum will stay as low is 
in Fig. El Such systematic variations may occur, for example, if sunspots and plages are 
not randomly distributed on the solar surface, but tend to appear in confined areas called 
"active longitudes" in the literature. Currently, there seems to be no strong evidence for the 
existence of such long-term structures. A possible north-south asymmetry in the distribution 
of magnetic features evolving with the solar cycle is another interesting route of investigation 
where the reconstructed TSI images can be used. 

The results derived in this paper apply to solar-type stars seen equator-on (inclination ~ 
90°). A pole-on configuration is optimal for astrometric detection, because both dimensions 
are equally engaged, but it is statistically less probable. Even at moderate inclinations to 
the line of sight, the distribution of magnetic perturbations in y-coordinate may be different 
from that we find for the Sun, if the magnetic features are confined to the equatorial zone. 
The features would normally be seen closer to the side of the limb opposite to the visible 
pole, producing a skewed distribution of A^ depending on whether the spots or bright areas 
are the dominating contributors. However, giant near-polar spots appear to be common on 
very active, fast-rotating stars. 

The research described in this paper was carried out at the Jet Propulsion Laboratory, 
California Institute of Technology, under a contract with the National Aeronautics and Space 
Administration. Observations at the 150-foot tower telescope on Mt. Wilson have been 
supported over the years by grants from the National Aeronautics and Space Administration, 
the National Science Foundation (NSF) and Office of Naval Research. Research is currently 
supported by the NSF through grant AGS-0958779. 
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